function [snear] = samplesAdjacency_fem (  x_samples , conectivity, gPts) %, order

sPts = size(x_samples,1);
nElements = size (conectivity,1);

snear = {zeros(sPts,1)};

ig=0;

for elem= 1:nElements
    e_near    = conectivity(elem,:); %vertices
    for gp =1 : gPts
        ig=ig+1;
        snear {ig} = e_near;
    end
end

end